Ultra-High-Density Molecular Core and Warped Nuclear Disk in the 
Deep Potential of Radio-Lobe Galaxy NGC 3079 
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ABSTRACT 

We have performed high-resolution synthesis observations of the 12 CO (J=l-0) line 
emission from the radio lobe edge-on spiral NGC 3079 using a 7-element mm-wave 
interferometer at Nobeyama, which consisted of the 45-m telescope and 6-element Array. 
The molecular nuclear disk (NMD) of 750 pc radius is found to be inclined by 20° from 
the optical disk, and the NMD has spiral arms. An ultra-high-density molecular core 
(UHC) was found at the nucleus. The gaseous mass of the UHC within 125 pc radius 
& ■ is as large as ~ 3 x 10 8 M®, an order of magnitude more massive than that in the same 

area of the Galactic Center, and the mean density is as high as ~ 3 x 10 3 H2 cm -3 . A 
position- velocity diagram along the major axis indicates that the rotation curve starts at 
a finite velocity exceeding 300 km s _1 already from the nucleus. The surface mass density 
in the central region is estimated to be as high as ~ 10 5 M(g pc -2 , producing a very deep 
. gravitational potential. We argue that the very large differential rotation in such a deep 

potential will keep the UHC gravitationally stable from current star formation. 

Subject headings: ISM: molecules — galaxies: individual (NGC 3079) — galaxies: 
ISM — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: structure - 
radio lines: ISM 
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1. Introduction 

NGC3079 is an edge-on galaxy known for its pro- 
nounced radio, Ha and X-ray lobes emerging from 
the nucleus (Hummel et al 1983; Duric et al 1983; 
Lord et al 1986; Veilleux et al 1994; Pietsch et al 

1999) . The nuclear outflow may originate with ac- 
cretion of dense gas onto a compact core (Irwin and 
Seaquist 1988), and the nucleus exhibits LINER and 
Seyfert 2 nuclear activity (Lord et al 1986). NGC 
3079 has been studied in various radio wavebands in 
continuum, HI and CO lines as well as maser lines 
(Irwin and Seaquist 1991; Ford et al. 1986; Duric et 
al 1983; Duric and Seaquist 1988; Irwin et al 1987; 
Irwin et al 1988; Young et al. 1988; Sofue and Ir- 
win 1992; Irwin and Sofue 1992; Sawada-Sato et al. 

2000) . The galaxy's distance is 15.6 Mpc (H = 75 
km s _1 /Mpc), and the major axis position angle of 
the optical disk is 167° and the inclination angle 85°. 

Earlier CO observations with the Nobeyama mm 
Array (NMA) inferred the existence of a dense nuclear 
molecular disk, rotating at a high velocity within the 
central 10" (750 pc) (Sofue and Irwin 1992; Irwin 
and Sofue 1992). NGC3079 is a galaxy whose un- 
usual nuclear activity can be probed by observing the 
affected interstellar medium, and can be compared 
with the Milky Way for its similar edge-on orienta- 
tion. Knowledge of the kinematics, morphology and 
density of the interstellar medium in the immediate 
vicinity of nucleus is particularly important, since the 
ISM can both fuel the nucleus and trigger activity. 
The nuclear disk would also play a role in collimating 
outflows. For this, higher resolution CO observations 
have been required to determine how the morphol- 
ogy and kinematics of the nuclear molecular disk are 
related to the nuclear activity. 

In this paper we report the result of high-resolution 
CO-line observations of NGC 3079 using a 7-element 
millimeter wave interferometer at Nobeyama by com- 
bining the NMA (Nobeyama mm Array) and 45-m 
telescope, whose code name was "RAINBOW". 

2. Observations 

The 12 CO (J=l-0) observations of NGC3079 were 
made on 2000 January 18 and February 2 in the 
7-element "RAINBOW" mode at Nobeyama, which 
consisted of six 10-m antennae array in AB configu- 
ration linked with the 45-m telescope. Observations 
in the C and D-configuration observations were made 
on 2000 March 12 and April 12, respectively. The 



UV coverage of the observations was ideal for the 
high declination. After obtaining UV data for indi- 
vidual array configurations, they were combined, and 
data analysis was done by using the AIPS standard 
packages. The UV data were CLEANed and Fourier 
transformed to three dimensional cubes in RA, Decl, 
and frequency space, and further transposed to (RA, 
Decl, Vlsr ) cubes. The synthesized HPBW was 
1".62 x 1".34, slightly elongated in the direction of 
position angle 110°, for combined UV data of the 
7-element synthesis in AB, C and D configurations 
(RABCD configuration). The beam for the 7-element 
AB configuration without C and D was 1".27 x 1".09 
elongated at a position angle 106° (RAB configura- 
tion). At a distance of 15.6 Mpc 1" corresponds to 75 
pc. 

The center position of the galaxy was taken at the 
radio nucleus at RA = 09h 58m 35.02s and Dec = 
55°55'15.4" (epoch 1950) and the systemic LSR ve- 
locity was taken as 1 1 13 km s~ 1 . The phase and band- 
pass calibrations were made by observing the nearby 
radio source QSO 0957+561, which was measured to 
have a flux density of 0.68 Jy at the observing fre- 
quency for RAB array in January and February, and 
0.58 Jy for C and D arrays in March and April 2000. 
The field of view for the 7-elelmcnt array including 
the 45-m telescope (RABCD) was 20", and that for 
the 6-element ABCD-configuration array without 45- 
m telescope was 60". This yielded higher sensitivity 
in the central 20" region than in the surrounding re- 
gion. 

Fig. 1 shows the integrated intensity map of the 
CO line from the RABCD configuration, and a mean- 
velocity field at a resolution of 1".62 x 1".36 is shown 
as the inset. The major structures discussed below 
are indicated by the arrows and dashed illustrations. 
A DSS S-band image is inset in the upper-left corner. 
Fig. 2 shows an integrated intensity for the central 
region from the RAB configuration with a higher reso- 
lution of 1". 27 x 1". 09. Fig. 3 shows position-velocity 
(PV) diagrams along the major axis with a 5" slit 
width by RABCD array (upper panel), and for the 
innermost region with 2" slit width by RAB array 
(lower panel). 

— Fig. 1, 2, 3 — 

3. Molecular Structures 
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3.1. Ultra-High-Density Molecular Core 

The most remarkable feature in the intensity maps 
(Fig. 1 and 2) is a very compact, intense CO concen- 
tration at the nucleus. This "ultra-high-density core 
(UHC)" of molecular gas is elongated in the north- 
south direction at PA=176°, inclined by 9° from the 
main disk. The full width of half maximum of the 
UHC is measured to be 3".3 x 1".8 (250 x 180 pc), 
and therefore, the radius is about 125 pc. 

The integrated CO intensity at the nucleus is 117 
Jy/beam km s _1 , or I C o = 5.0 x 10 3 K km s _1 . The 
H 2 column density toward the nucleus is estimated 
to be Nn 2 ~ 5.0 x 10 23 H 2 cm~ 2 , where we took a 
conversion factor for the centers of galaxies of X = 
1.0 x 10 20 H 2 K-^km s" 1 ]" 1 (Arimoto et al 1996). 

The total molecular mass in the UHC is estimated 
to be M gas = 3x 10 8 M m . Since the vertical direction 
is not resolved, we here assume that the thickness of 
the UHC is of the same order of 30 pc as the molecular 
ring in the Milky Way center (Sofue 1995). Then, the 
mean density of molecular hydrogen is estimated to 
be of the order of or higher than ~ 3 x 10 3 H 2 cm~ 3 . 
Hence, the UHC consists of a pile of molecular gas 
equivalent to a thousand medium-sized GMCs within 
a 125 pc radius. 

The velocity field in Fig. 1 (lower-right inset) 
shows that the UHC is rotating regularly but very 
rapidly. The PV diagrams (Fig. 3) show that the 
rotation velocity in the central 1" region rises ex- 
tremely steeply, or more likely, the velocity starts at 
a finite velocity of ~ 300 km s _1 already at the nu- 
cleus. We may estimate the dynamical mass from 
the rotation velocity because the pressure term due 
to velocity dispersion is negligible for population I 
gases, particularly for molecular gas. The dynami- 
cal mass within R = 125 pc is, then, estimated as 
M dyn = V r 2 ot R/G - 2 x 10 9 M(g. Hence, the gaseous 
mass in the molecular core shares about 15% of the 
dynamical mass. 

3.2. Warped Nuclear Disk with Spiral Arms 

The molecular core is surrounded by the nuclear 
molecular disk (NMD) of radius 10" (750 pc) elon- 
gated in the direction of the optical major axis (dashed 
ellipse in Fig. 1). The projected minor-to-major axial 
ratio of the ellipse is as large as 0.40. Since the disk 
thickness is considered to be sufficiently small, such 
a fat shape indicates that the inclination angle of the 
disk is either 66°. This implies that the NMD is not 



in the same plane as the outer optical disk whose in- 
clination angle is 85°. The NMD is, therefore, warped 
by about 20° from the main optical disk. 

The disk is superposed by two symmetrical spiral 
arms, one arm extending from the north-west of the 
UHC toward the north, and the other from south-east 
to the south, which are illustrated by dashed lines in 
the inset in Fig. 1 . If we assume that the spiral arms 
are trailing, the eastern half of the disk is in the far 
side, and western half is near side. This configuration 
is in the same sense as that of the outer optical disk. 

The velocity field shows that the circular rotation 
is dominant. The slightly distorted velocity field is 
superposed, which suggests small non-circular motion 
due to spiral arms. Since the intensity map shows 
spiral arms, the latter may be more likely In the PV 
diagram the spiral arms show up as inclined ridges in 
the central 10", crossing the bright ridge of the UHC. 

The total mass of molecular gas of the NMD in- 
cluding the outskirt extending to radius 900 pc is es- 
timated to be 2 x 10 9 M ( g ) for the same conversion fac- 
tor as above. This shares about 14% of the dynamical 
mass within the same radius, 1.4 x 10 10 M,g, for a ro- 
tation velocity of 260 km s _1 at 900 pc radius. 

3.3. Main Disk 

Fig. 1 indicates the existence of a more extended, 
narrower component extended in the direction of the 
optical major axis for more than ±30" (±2.2 kpc). 
This feature is attributed to the outer main disk visi- 
ble in optical images associated with dust lanes. This 
disk component is slightly lopsided in the sense that 
the southern part is brighter. 

The disk shows up also in the PV diagram (Fig. 
3) as the two symmetrical tilted ridges crossing the 
UHC and NMD. Such PV ridges are typical for grand- 
designed two armed density-wave spirals. These PV 
ridges are bifurcated from the maximum- velocity parts 
of the PV ridges due to the NMD. Hence, these disk 
arms may be continued to the spiral arms in the 
NMD. 

4. Rotation Curve and Mass Distribution 

Applying the envelope-tracing method (Sofue 1996) 
to the PV diagrams, we have derived a rotation curve, 
and combined it with an outer rotation curve obtained 
from HI data (Irwin and Seaquist 1991). The ob- 
tained rotation curve for the entire galaxy is shown 
in Fig. 4. The rotation curve starts already from 
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a very high finite value in the nucleus, followed by 
a usual curve in the disk and halo. High velocities 
in the center, followed by a flat rotation in the disk 
and halo, exhibiting a broad maximum in the disk, 
are commonly found in spiral galaxies observed with 
sufficiently high resolutions (Sofue et al. f999). 

Since the circular motion is dominant (Fig. 1), we, 
then, use the rotation curve to directly calculate the 
distribution of surface-mass density (SMD), applying 
the method developed by Takamiya and Sofue (1999). 
Fig. 5 shows the obtained SMD, where the thick line 
shows the result for a flat-disk assumption, and the 
thin line for a spherical assumption [see Takamiya and 
Sofue (1999) for detaileds] . The true SMD value is 
considered to be sandwitched by the two lines (thick 
and thin) within an error of about 30%, except for 
the outermost region, where the edge effect becomes 
not negligible. The SMD between R ~ 1 and 8 kpc is 
well fitted by an exponential disk of a scale (e-folding) 
radius of 3.5 kpc. The bulge component between 0.3 
and 1 kpc can be approximated by an ~ e -( R / 20 °P c ) 1 
law. The innermost region within R ~ 300 pc shows a 
much steeper increase toward the nucleus, exhibiting 
an extremely high-density dynamical core with SMD 
as high as l0 5 M m pc~ 2 at R ~ 100 pc. 

— Fig. 4, 5 - 

5. Discussion 

The central R ~ 100 pc region of NGC 3079 is 
filled with the UHC with a molecular mass as high as 
3 x 10 8 M®. The UHC is rotating regularly at high 
speed of ~ 300 km s _1 in the deep gravitational po- 
tential. The molecular gas density and mass are an 
order of magnitude greater than those in the Galactic 
Center. We stress that such an extremely high den- 
sity, massive interstellar condition was found for the 
first time. 

We here consider the reason why such a dense core 
could have survived in the gas phase, without suffer- 
ing from current star formation. The Jeans time in 
the UHC is of the order of tj ~ 10 6 yr. On the other 
hand, the dynamical time scale for a cloud to be torn 
off by the Coriolis force and differential rotation is 
t u ~ dV/dR - w~2.5x 10 5 yr. Here, V, R and w 
are the rotation velocity, radius, and angular veloc- 
ity. Hence, tj > tu, so that clouds cannot collapse to 
form stars, being are kept gravitationally stable and 
stretched azimuthally along the orbits by the differen- 
tial rotation. Velocity dispersion would be an alterna- 



tive cause for surpression of star formation. However, 
the random motion and turbulence in molecular gas 
are known to be several km s _1 or less, which will less 
affect the stability compared to the effect of the high 
differential rotation. 

Formation of the UHC will be understood as due to 
accretion by non-circular streaming motion from the 
NMD, as readily suggested from the spiral pattern in 
the NMD. The existence of a compact, massive central 
mass would cause rapid accretion of the disk gas to a 
dense core in circular rotation, as is indeed simulated 
by Fukuda et al. (2000). 

Another important finding is that both the NMD 
and UHC are significantly warped from the main disk. 
The warp of NMD may explain why the Ha lobe is 
highly asymmetric, or lopsided (Veilleux et al 1994). 
Warping of the central gas disks appears not a rare 
case. In fact, Schinnerer et al. (2000) have shown 
by their high-resolution interferometer observations 
in the CO line, that the nuclear molecular disk of the 
Seyfert galaxy NGC 1068 is significantly warping. We 
finlly mention that the origin of the highly inclined 
nuclear torus of a few pc scale, as inferred from VLBI 
jet directions (Irwin and Seaquist 1988) and maser 
spots (Sawada-Satoh et al. 2000), would be somehow 
related to the warp of its progenitor, NMD and UHC, 
from which the torus had accreted. 

The authors thank the NRO staff for helping us 
in the observations. They also thank T. Takamiya 
for the program to calculate the surface mass density 
from rotation curves. 
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Figure Captions 

Fig. 1. Distribution of integrated intensity of 
the 12 CO (J = 1 — 0) line emission in the central 
40" x 40" area of NGC 3079. The angular resolu- 
tion is 1".62 x 1".34. An ultra-high-density molecu- 
lar core (UHC) is embedded in the nuclear molecu- 
lar disk (NMD), which are further surrounded by the 
main disk component (arrows). Contours are at 2.5, 
5, ...17.5, 20, 25, 35, 40, 50, ... 90 and 100% of the 
peak intensity. 

Lower-right inset is a mean-velocity field. The 
thick velocity contour near the center is at 1100 
km s _1 , and the contours are drawn at 50 km s _1 
inclement, increasing toward the south. The thin 
dashed ellipse outlines the NMD. Two molecular spi- 
ral arms are traced by the thick dashed lines. 

Upper-left inset is a £?-band image from DSS (ver- 
tical extent 8'). The small white box indicates the 
CO map region. 

Fig. 2. The central 20" x 20" at a higher resolution 
of 1".27 x 1".09. The UHC shows up more clearly. 
Contours are at 5, 10, 20, ... 90, and 100% of the 
peak intensity. 

Fig. 3. The upper panel shows a position- velocity 
diagram along the major axis of NGC 3079 in CO 
with a slit width of 5 arcsec (original resolution 
1".62 x 1".34). Contours are at every 5, 10,..., 25, 
30, 40, .., 90 and 100% of the peak intensity. The 
lower panel shows the same, but for the central re- 
gion at higher resolution with a slit width 2 arcsec 
(original resolution 1".27 x 1".09). 

Fig. 4. The rotation curve of NGC 3079. 

Fig. 5. Surface mass density plotted against ra- 
dius. The thick curve is the result for a flat-disk as- 
sumption, the thin line for a spherical assumption. 
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